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New rotation laws have been recently found for general-relativistic self-gravitating stationary 
fluids. It was not clear whether they apply to systems rotating with a constant linear velocity. In 
this paper we fill this gap. The answer is positive. That means, in particular, that these systems 
should exhibit the recently discovered general-relativistic weak-field effects within rotating tori: the 
dynamic anti-dragging and the deviation from the Keplerian motion induced by the fluid selfgravity. 


INTRODUCTION 


Axially symmetric and stationary Newtonian hydro- 
dynamic configurations are known for a long time to be 
characterized by a rich variety of rotation curves. The 
angular momentum per unit mass j can be any func¬ 
tion of r, where r is the distance from a (fixed) rotation 
axis. In contrast to that, in general relativity had been 
known only two families of rotation laws—one with 7 
being a linear function of the angular velocity 17 [iHsI 
and a more recent nonlinear angular velocity proposal 
Q. Their Newtonian limits recover only a small fraction 
of the set of Newtonian rotation curves. 

Quite recently two of us have found general-relativistic 
rotation curves j = @ that in the nonrelativistic 

limit exactly coincide with all monomial rotation laws 
ilo = with the exception of the constant lin¬ 

ear velocity case (—00 < <5 < 0, J 7 ^ —1, w = const). 
We obtained in particular the general-relativistic Keple¬ 
rian rotation law that possesses the first post-Newtonian 
limit (IPN) and exactly encompasses the solution corre¬ 
sponding to the massless disk of dust in the Schwarzschild 
spacetime. 

The rotation law proposed in Q reads 


jm = 




1 - 


1 - 3(5 

(U-(5)c2 





( 1 ) 


The main purpose of this paper is to show that the 
problematic case of constant linear velocity is also in¬ 
cluded in the proposed family of general-relativistic ro¬ 
tations. 


HYDRODYNAMICAL EQUATIONS 

In this paper we apply the formulation of general- 
relativistic hydrodynamics elaborated by Komatsu et al. 
in Q. Einstein equations read 

~ ( 2 ) 

where is the stress-energy tensor. We shall assume 
axial symmetry, stationary rotation and the angular ve¬ 
locity vector field v = Vld^. Then the metric can be 


written as 

ds^ = —e~^ {dx^Y + (^(j) — (r, z) dx°^ 

+e^ {dr^ +dz^), 

where r, z, (p are the cylindrical coordinates. The stress- 
energy tensor of a relativistic perfect fluid reads 

= p{c^ + h)u^u^ + . 

Here p is the baryonic rest-mass density, h is the spe¬ 
cific entalpy and p is the pressure. The 4-velocity is 
normalized, gapu°‘u^ = —1 and /u*' = S7. 

We assume that the equation of state obeys a poly¬ 
tropic relation 

p{p,S)=K{S)p\ 


where S is the specific entropy of the fluid and 7 is known 
as the adiabatic index. Hence h{p,S) = K{S)^:^p'^~^. 
The entropy is constant. 

We introduce, following Q, 

where V is the proper velocity with respect to the zero 
angular momentum observer. 

Einstein equations implied by the above metric take 
the form of an overdetermined, but consistent set of equa¬ 
tions imposed on the potentials a, (3, v, and uj. These 
equations have been found by Komatsu et al. in Q. 

If we assume that the angular momentum per unit 
mass. 




J = Ufj,U = 




( 3 ) 


depends only on the angular velocity Q [j = j (S7)], then 
the Euler equations become solvable and they reduce 
to a single general-relativistic integro-algebraic Bernoulli 
equation 

In (l + ^)+^ + ^ln(l-^)+^ j dnm)=c, 

( 4 ) 

where C is an integration constant. The above equation 
carries all information that is present within the con¬ 
servation equations = 0 and the baryonic mass 

conservation [pu^) = 0 . 
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ROTATION LAW 


The general-relativistic rotation law employed in Q 
equates the angular momentum per unit mass, given by 
to a specific function 


jm = 


w 




1 1-3(5 „,,i-,5r)l-l-(5 I 4co 

i - JT+SJ^W 


( 5 ) 


In explicit terms one has 




y 2 




( 6 ) 


From this equation one can recover rotation curves 
n{r,z). 

With the rotation law m the general-relativistic 
Bernoulli equation (|T|) acquires a simple algebraic form, 
assuming S ^ —1: 


\ ( 1 - 35 ) 

j = C. (7) 

Assume that there exists the Newtonian limit (the ze¬ 
roth order of the post-Newtonian expansion—OPN) of 
the rotation law. This yields, in the OPN order. 


1 + F 1 e-/'- 


1 - 


1-3(5 
;2(l + ,5) 




Define uj = r ‘^A^. The spatial part of the metric 
ds^ = — -I- ^ -I- {dx°y — 2c~^A^dx^d(f) 

{dr^ + dz^ + r^d(j?) . (9) 

is conformally flat. 

We split different quantities (p, p, h, U, and u®) into 
their Newtonian (denoted by subscript ‘0’ and IPN (de¬ 
noted by subscript ‘1’) parts. Exempli gratia, for p, D, 
'll, and U this splitting reads 



p = po + c~‘^pi, 

( 10 a) 


U = Ug + , 

( 10 b) 


u = uo + c-y. 

( 10 c) 

Notice that, 

, up to the IPN order. 



-dip = diho + c~‘^dihi+0{c~y, 

P 

( 11 ) 

where the IPN correction hi to the specific enthalpy can 
be written as hi = j^Pi- For the polytropic equation of 
state this gives hi = (7 — 1 ) hopi/po- 

Making use of the introduced above splitting of quanti- 


ties into Newtonian OPN and IPN parts one can extract 
from Eq. © the OPN- and IPN-level Bernoulli equations. 

At the Newtonian level the gravitational potential is 
given by the Poisson equation 


1 -^ 


^^0 = ( 8 ) 
j' 1 — 5 

Thus w and S can be obtained from the Newtonian limit. 
Moreover, the constant w is any real number, while S is 
nonpositive—due to the stability requirement Q— and 
satisfies the bounds —oo < 5 < 0 and 6^—1. These two 
constants can be given apriori within the given range 
of values. Let us point out that the general-relativistic 
extension of this stability condition, formulated in 0 , is 
satisfied in our case for 5 < 0. 

The two limiting cases 5 = 0 and 5 = — oo correspond 
to the constant angular momentum per unit mass (fig = 
w/r'^) and the rigid rotation (f2 = w), respectively. The 
Keplerian rotation is related to the choice of 5 = —1/3 
and = GM, where M is a mass. 

The particular form of the expression (1 — 35)/(l -I- 5) 
in the denominator of ® follows from the condition that 
a rotating infinitely thin disk made of weightless dust in 
a Schwarzschild space-time satisfies exactly the Bernoulli 
equation and the Kmlerian rotation law [^. 

It was proven in y that if cg is the Newtonian hydro- 
dynamic energy per unit mass, then the exact solution 
satisfies the first post-Newtonian (IPN) equations. We 
shall outline here the calculation. 

The IPN approximation corresponds to the choice of 
metric exponents a = (3 = —v = —U, with \U\ <C Q. 


AUg = 4^Gpg, 


while the Bernoulli equation reads 


ho + Uo 


5-1 
2(1 + 5) 


ny 


= Co, 


( 12 ) 


(13) 


where cg is a constant that can be interpreted as the 
energy per unit mass. Here A is the flat Laplacian with 
respect to the cylindrical coordinates r, z, and 4>. 

One can obtain the first correction vf to the angular 
velocity H by expanding the rotation law (jH]) in powers 
of c up to terms (c“^) 


= - 


1-5 


ny + 




0/^0 


r2(l-5) 1-5 


(14) 


Using the fact that in the Newtonian gauge imposed on 
the line element (O the distance to the rotation axis mea¬ 
sures f = r (l — Uo/y + 0{c~y, we can write down 
the full expression for the angular velocity, up to terms 

oic-y i 


4> 

n = H—T — 

UJ 2 

f2/(l-S) ~ c2(l _ 5) ^0 (^0 + ^0^^) + 

A(j) 4 

+ f 2 c 2 (l_J) “ c2(l-5)^°^°' 


( 15 ) 
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CONSTANT LINEAR VELOCITY 


In order to construct the limit (5^—1, we need to in¬ 
spect the denominator of One can easily check, using 
([TS]), that in this limit both cq and e2^(7+g) be¬ 

come singular. After addition and subtraction of a term 
c‘^{i+s) the denominator of ([S]) we obtain 


= 




1 1-35 „„l-5r)l-l-5 2(5-1) , 4co 

~ (l+ 5 )e^ ^ ^ + T3- 


^(T+5y' 


W 




1 + y„.-<ni+i _ 


( 16 ) 


where cq is the regularised energy per unit mass 


co = ho + Uo- A 


(5-1 


,.1-5 


— Hq + Uq — 


2(1 + 5) 2(1 + 5) 

5-1 

V 1 + 5 


-w 


(17) 


It is worth noting that such defined cq and j(i2) are not 
singular in 5 = — 1. To prove this fact, one needs only to 
use the identity 

x“ — 1 

lim- = In a;. 

Q:—>-0 CH 

Indeed, taking the limit <5 —>■ — 1 in equations m and 
(HZl) results in 


m 


2o-l 




l + 4w^2_ 4^2 


4co 


while the regularised Bernoulli equation reads now 

Co =/iQ + In n. (18) 

As in the case 5 — 1, the leading correction vf to the 

angular velocity Hq is obtained from the perturbation 
expansion of the rotation law 




W ^ In O _j_ 4cq 

^2 ^2 I (.2 


V 2 

(1 


C 2 ) 


(19) 


the geometric distance, as given at the IPN level of ap¬ 
proximation, from the rotation axis. The second term 
in (HID is sensitive both to the contribution of the disk 
self-gravity at the plane z = 0 and the deviation from the 
strictly Keplerian motion. It vanishes for test fluids. The 
third term is responsible for the geometric frame drag¬ 
ging. The last term represents the recently discovered 
dynamic anti-dragging effect; it agrees (for the monomial 
angular velocities flo = —with the result ob¬ 

tained earlier in Q. 

We shall consider the remaining first order perturba¬ 
tions terms. 

The vectorial component satisfies the following 
equation 


AA^ - 2 


dr A, 


— = —WnGr'^poilo. 


( 21 ) 


The IPN Bernoulli equation has the form 


Cl = -hi -Ui- floA^ + 2r^(no)^^o - 


-4hoUo - 2t/2 - 


4(1 + 5) 


— —hi — Ui — ^IqA^ + 2r^(f2o)^^o ~ 2^0 

4 


-4/ioC/o - 2U^ - 


(5 + 3)w' 
4(1 + 5) ’ 


( 22 ) 


where ci is a constant. 

The IPN potential correction Ui can be obtained from 

AUi = AttG {pi + 2po + poiho - 2Uo + 2r^ino)^)) . 

(23) 

It is clear that only the correction to the proper en¬ 
ergy per unit mass ci would cause trouble in the limit 
5 —>■ — 1, but on the other hand it does not influence the 
IPN correction to the angular velocity. One would have 
to regularise ci in the next orders of the perturbation 
calculation. 


SUMMARY 


up to terms of the order c ^. One arrives at 

vt = -ny + ^ - 2O0/10, ( 20 ) 

where we applied Eqs. dS]) and m- Note that the form 
of the corection is not affected by the transition to the 
limit and agrees with (fHl) . 

After these consideration we are able to interpret the 
meaning of various contributions to the angular veloc¬ 
ity fl given by formula (1151) . The first term is simply 
the Newtonian rotation law rewritten as a function of 


We found in Q that the consistency of the formalism 
and the uniqueness of those solutions that possess the 
1st post-Newtonian limit, lead to well defined numerical 
values of the coefficients in the rotation law O- It is in¬ 
teresting that the definition is rigid—there is no any 
parametric freedom left for those solutions that have the 
1st post-Newtonian limit. We proved in this paper that 
the case of constant linear velocity, that corresponds to 
the seemingly singular point 5 = — 1, can be deduced 
from ©• That extends the validity of the general rela¬ 
tivistic laws constructed in Q. 
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